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ABSTRACT 

We have made a sensitive search for the origin bands in the known electronic transi- 
tions of the linear carbon chains C4 and C5 at 3789 and 5109 A towards £ Oph (Ay<l). 
The incentive was a recent detection of C3 in this interstellar cloud with a column den- 
sity of 1.6xl0 12 cm -2 plus the availability of laboratory gas phase spectra of C4 and 
C5. Further, some models of diffuse interstellar clouds predict that the abundance of 
these latter species should be within an order of magnitude of C3. Despite achieving 
S/N of 2300 to 2600 per pixel at a resolution of ~110,000, the searches were negative, 
leading to 3<t upper limits to the column density of N(Cs)= 2xlO n cm -2 and N(C4)= 
4xl0 12 ~ 13 cm~ 2 where these values rely on theoretically calculated oscillator strengths. 
The implication of these limits are discussed on the choice of molecules for study in 
future attempts to identify the carriers of the stronger diffuse interstellar bands. 

Subject headings: ISM: molecules— C4, C5 



1. Introduction 

Carbon chain molecules are often at the forefront in discussions of the diffuse interstellar bands 
(DIB) (see Douglas (1977) and Smith et al. (1977)), which are found mainly in the optical part of 
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the spectrum (Herbig 1995). They became more appealing candidates with the discovery at mm 
wavelengths (McCarthy & Thaddeus 2001) of numerous polar carbon chains in dense interstellar 
clouds. However, only in the past few years have the gas phase electronic spectra of a number of 
such species (e.g. C4, C5, C2 n H n=3-6), as well of related cations (e.g. HC2riH + n=2-4, HC2 n N + 
n=2-4) and anions (e.g. C~ n=3-ll) been detected in the laboratory, enabling direct comparison 
with astronomical data (for example, Motylewski et al. (1999)). In all cases where the comparisons 
could be made, the results were negative, leading to the conclusion that the column densities of 
these species are < 10 12 cm~ 2 in diffuse clouds. 

However, with the detection of the rotational lines in the electronic transition of C3 near 4052 
A (Maier et al. 2001), corresponding to total column densities of l-2xl0 12 cm~ 2 , and the prediction 
of certain models of such diffuse clouds (e.g. Terzieva & Herbst (1998)) that the abundances of 
C4 and C5 are less than a factor of ten smaller than that of C3, it became appealing to search 
for these molecules. A disadvantage compared to the electronic spectrum of C3 is that, according 
to theoretical predictions, the oscillator strengths of the transitions appear to be smaller. Unfor- 
tunately, experimentally determined oscillator strengths are unavailable. These smaller oscillator 
strengths result in higher abundance limits than was the case for C3. Apart from the abundances 
predicted for C4 and C5, the detection of both C3 and C5 in the infrared by Bernath et al. (1989) 
in a circumstellar shell provides a guideline. In the latter work the concentration of C5 proved to 
be some ten times less than that of C3. 

Here we report an attempt to detect the origin bands of the electronic transitions of C4 ( 3 S~ 
— 3 £~ ) at 3789 A and of C5 ( 1 n u — ) at 5109 A in absorption towards the reddened star 
C Oph (HD 149757) with the Gecko spectrograph of the Canada-France-Hawaii 3.6 m telescope. 
There is already one report in the literature of a non-detection of C5 (Galazutdinov et al. 2001). 
The present study is an order of magnitude more sensitive. There is no prior attempt reported to 
detect C4 in interstellar clouds, but the electronic spectrum in the gas phase was only obtained a 
year ago (Linnartz et al. 2000). 

2. The Observations 

The reddened star ( Oph (HD 149757), was observed on 29 and 30 June 2001 (UT) with the 
Gecko echellette spectrograph, fiber fed from the Cassegrain focus of the Canada-France-Hawaii 3.6- 
m telescope (CFHT) (Baudrand & Vitry 2000). This star is bright, having a visual extinction, A„, 
near 1 with a rich spectrum of sharp interstellar lines. Crawford (1997) has resolved the interstellar 
C2 at 8756 A into two close velocity components separated by 1.1 km s _1 which is hard to resolve 
at our resolution. The rapidly rotating star, rj UMa (HD 120315) was observed as standard. 

The detector was the rear illuminated EEV1 CCD (13.5 //m 2 pixels) and the spectral regions 
were centered at 3789 A in the 15th order, and at 5109 A in the 11th order. The 15th order was 
isolated by the Gecko ultraviolet grism, the 11th by the blue grism. Individual spectra of ( Oph 
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were exposed for 15 minutes at 3789 A and 6 minutes at 5109 A. The feed fibre was continuously 
agitated to overcome modal noise (see Baudrand & Walker (2000)). Lines in Th/Ar comparison 
arc spectra, taken before and after each set of stellar spectra, typically had FWHM of 3 pixels, 
corresponding to resolutions of R = 115000 and 109800 at 3789 and 5109 A, or 0.0109 and 0.0155 
A pxl -1 , respectively. Extensive series of flat-field spectra of a quartz-iodide lamp were recorded 
for each spectrograph setting and groups of biases were taken several times each night. 

Conditions were not ideal to achieve the very high S/N (>4000) which we felt necessary to 
detect the C4 and C5 bands. Only parts of two of the assigned nights were clear and the flat-field 
spectra displayed an unstable structure typical of modal noise. The stellar spectra did not show 
the same structure and we were able to demonstrate that the main fiber agitation was working as 
designed but a misalignment of the fiber feeding the the light of the flat field lamp to the main 
fiber may have caused the problem. For this reason a combination of flat fields and standard stellar 
spectra were used in the reduction. 

The many biases taken throughout the observing run were averaged to remove the zero-level 
offset in the spectra for both wavelength regions. For the 5109 A spectra, obtained after the noise 
problem with the flat field exposures had been recognized, rj UMa was observed at a S/N level 
substantially higher than £ Oph. Since the spectrum of the standard star is featureless near 5109 
A its spectrum was used in place of flat field spectra. 

Unfortunately, at 3789 A, weather constraints prevented observations of 77 UMa at a sufficient 
S/N level for them to be used directly as flat fields. Instead, we averaged the flat field exposure 
to remove pixel-to-pixel sensitivity variations in the spectra of both i] UMa and £ Oph. One 
dimensional spectra of both stars were then extracted in a standard manner. Because of the modal 
noise in the flat field spectra, each flat field corrected stellar spectrum shows a residual low frequency 
modal noise pattern. We therefore smoothed the featureless, flat fielded spectrum of 77 UMa with 
a 10-point box car filter and then divided the ( Oph spectrum by this smoothed spectrum to 
remove the modal noise residuals without seriously degrading the inherently high S/N of the £ Oph 
spectrum. For clarity this reduction procedure is illustrated in Figure 1. 

Low-order polynomial fits to the positions of the Th/Ar arc spectra were used to calibrate the 
£ Oph spectra in wavelength. The spectra were then normalized to the continuum and heliocentric 
corrections applied. The observations are summarised in Table 1 which lists exposure times and 
S/N per pixel for each spectral region. The final column gives the radial velocities of the interstellar 
K I 4044.1 and 4047.2 A lines as quoted by Maier et al. (2001). These velocities were applied to 
each spectrum to put the interstellar features on a laboratory scale before making the comparisons 
discussed in the next section. The comparisons are shown in Figures 2 and 3 where the stellar 
spectra have been smoothed with a 3-pixel box car filter. 

The 3<7 detection limits are derived from: 



W max = 3(wd) 1 HS/N)- 1 
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where the 3a limiting equivalent width, W max , and the FWHM of the feature, w, are both measured 
in A, the spectrograph dispersion, d, in A pixel -1 , and S/N is the signal to noise per pixel. We 
adopted w = 0.24 and 0.13 for the 3789 and 5109 A laboratory features, respectively. 

3. Results and Discussion 

Our attempt to detect the C4 3789 A origin band of the known electronic systems in absorption 
through the diffuse cloud towards the reddened star ( Oph was negative. In Figure 2 we compare 
the stellar spectrum with one from the laboratory (Linnartz et al. 2000) recorded at a temperature 
of around 50 K. This temperature should be representative for a non-polar molecule in diffuse 
interstellar clouds; in the case of C2 (Lambert et al. 1995) and C3 (Maier et al. 2001), comparable 
temperatures have been inferred. The laboratory spectrum (which differs from the published one 
only in that the overlapping C2 lines have been removed) shows that the rotational lines are lifetime 
broadened allowing only the lines in the P-branch to be resolved. The intense part of the band to 
the blue is the R-head and this feature was primarily sought in the astronomical spectrum. The 3a 
upper limit for the equivalent width (W max ) detection of the R-head is 0.67xl0 -4 A. The upper 
limit to the column density N max can then be derived from: 

N max = 1.13 x W 20 W max /X 2 f cm" 2 

where / is the oscillator strength of the band at wavelength A in A. 

Though the 3 S~ — 3 £~ electronic transition of C4 has been observed in the gas phase, the 
oscillator strength / is not known experimentally. Two values from ab initio calculations are avail- 
able which, however, differ by an order of magnitude. The first study yielded /=0.003 (Pacchioni & 
Koutecky 1988) whereas the second gave f=0.0005 (Muhlhauser et al. 2000). The values are for the 
whole band system and thus the /o-o for the 3789 A band will be reduced by the Franck-Condon 
factor for this transition. Based on the intensity distribution of the vibrational bands observed in 
the absorption spectrum measured in a neon matrix (Freivogel et al. 1996), this may be a factor of 
5. Thus, taking / _ in the 0.0001 to 0.001 range leads to N(C 4 ) < 4x 10 12 to 10 13 cm" 2 . 

Two more recent models of diffuse interstellar clouds (of which the data have been made 
available to the authors) with characteristics comparable to the one towards £ Oph (van Dishoeck 
& Black 1986) yield N(C 3 )/ N(C 4 ) ratios of about 2 (Terzieva & Herbst 1998), and 13 (Turner 
2000). The detected total column density of C 3 is 1.6xl0 12 cm' 2 (Maier et al. 2001). Thus, in 
view of the uncertainty in the oscillator strength, it is not possible to assess the predictive value of 
the models. 

The situation is similar for C5. The band system with the 5109 A origin band shown in Figure 3 
was observed first in absorption in a 5 K neon matrix (Forney et al. 1996) and then in the gas phase 
(Motylewski et al. 1999). It was assigned as the l Ii u — 1 S+ electronic transition by analogy to the 
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comet band system of C3. The rotational structure is only partially resolved due to line widths of 
0.7 cm -1 (homogeneous broadening due to intramolecular processes). The oscillator strength /o-o 
= 0.02 was estimated by taking the experimentally known value for C3 (/o-o = 0.016) and scaling 
this up by the length of the molecule, as simple quantum models such as particle in a box predict. 

However, on the basis of recent theoretical calculations, it is proposed that the 5109 A band 
system is in fact a forbidden transition, to 1 A M and/or states (Hanrath & Peyerimhoff 2001). 
The oscillator strength could not be calculated because the vibronic effects which lead to its intensity 
would have to be correctly accounted for. 

There is another band system in the neon absorption spectrum with origin band at 4454 A 
(unpublished data from the laboratory at Basel), which would then be the l U u - X S+ electronic 
transition, with calculated oscillator strength of 0.03 (Hanrath & Peyerimhoff 2001); an earlier f 
value for this transition was 0.037 (Kolbuszewski 1995). A pragmatic approach to the estimation 
of the oscillator strength for the 5109 A band is to use the relative intensities of the absorption 
systems in the 5 K neon matrix. The 4454 A absorption band system appears about a factor of 5 
more intense than the 5109 A one, implying an f value for the system of around 0.006. The /o-o 
value for the 5109 A band will then be reduced further by its Franck-Condon factor, to yield /o-o 
around 0.001. The lower spectrum in Figure 3 is that recorded towards ( Oph. The signal-to-noise 
ratio (2600) is high, leading to a 3a detection limit of 5xl0 -5 A. As the laboratory spectrum 
corresponds to a temperature of around 50 K, about 10 rotational lines comprise the Q-band head. 
This means that our astronomical measurements had a detection limit for an individual rotational 
line of ~ 5 x 10~ 6 A. For a lower temperature of say, 10 K, there would still be some five lines 
unresolved within the band to give a 3a limit per line of ~ 10" 5 A. 

Our 3a detetection limit of ~ 5 x 10~ 5 A for the band, together with /o-o = 0.001, leads to 
N(C 5 ) <2xlO n cm" 2 . Galazutdinov et al. (2001) recently reported a value of N(C 5 ) < 10 11 cur 2 , 
however the authors were unaware of the spectroscopic problems associated with the oscillator 
strengths and used a value of /o-o of 0.02. Taking the presently suggested /o-o =0.001 increases 
their value to N(Cs) < 2xl0 12 cm -2 , which is an order of magnitude less sensitive than our results. 

3.1. Conclusion 

According to the results reported here, N(C4) < 4xl0 12 to 10 13 cm -2 , and N(Cs) < 2xlO n 
cm -2 in the diffuse cloud towards £ Oph, while N(Cs) = 1.6xl0 12 cm~ 2 from the earlier study 
(Maier et al. 2001). From these column densities we calculate abundances relative to the total 
column density of hydrogen (n(H) + 2n(H2)), 1.3xl0 21 cm~ 2 (from Table 2 of van Dishoeck & 
Black (1986)), of C 3 = 2xl0~ 9 , C 4 <4xl0" 9 to lO" 10 , and C 5 < 2xl0- 10 . The results of the 
two diffuse cloud models to which we have access yield abundances which are too high for C3 by 
about an order of mag nitude, 3xl0~ 8 (Terzieva & Herbst 1998), and 5xl0~ 8 (Turner 2000). The 
values from the former model correspond to quasi steady state after 10 5 years and the latter are 
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for Ay = 1 and n= 500 cm - chosen with depletions which provide a good fit for essentially all 
molecular species in translucent clouds and for the majority of species in diffuse clouds. Similarly, 
the predicted abundances are also too high for the longer chains with relative values of CsiC^Cs 
of 5:3:8 from Terzieva & Herbst (1998), although Turner (2000) gives 25:2:1. So it appears that 
photodissociation rates have been underestimated (or other depletion mechanisms such as electron 
attachment or absorption on grains need to be included). The upper limit we have obtained of 
0.1 to the C5:C3 ratio is consistent with the result of infrared detection of both these species in 
a circumstellar shell with a ratio of 0.09 (Bernath et al. 1989), suggesting that our measurements 
may have been close to actually detecting C5. 

The detection of polyatomic species containing carbon atoms in diffuse and translucent inter- 
stellar clouds indicate comparable column densities, not larger than about 10 12 to 10 13 cm -2 . This 
set comprises the detection of C3 in the optical region (Maier et al. 2001) as well as of the rotational 
spectra of the polar species, C2H, C3H2 (Lucas k, Liszt 2000) as examples. So far among poly- 
atomic species only H^" has column densities in the diffuse medium exceeding 10 14 cm -2 (McCall 
et al. 1998). The implication of this for the search of appropriate molecular systems which could 
correspond to the stronger, narrower DIBs, with typical equivalent widths of 0.1 A, is as follows. 
Assuming that the species would have oscillator strengths of electronic transitions (in the visible) 
in the 1-10 range, as could be the case for carbon chains with 10-20 atoms, then the column density 
would have to be about 10 11 to 10 12 cm" 2 . This column density could be easily attained by the 
larger carbon species, whatever their shape, because they are less efficiently photodissociated. 

The suggestion by Douglas (1977) that long carbon chain molecules, C n (n = 5-15) be con- 
sidered as the carriers of DIBs, has now been tested directly for C4 and C5, showing that their 
abundance is too small. Thus the next step in attempting to identify carriers of the DIBs would 
be to obtain laboratory spectra of carbon species with 10-20 atoms having electronic transitions in 
the visible region (e.g., as is the case for C2n+i n=8-20, Maier (1998)) and then making a direct 
comparison with astronomical observations at high S/N as has been demonstrated in this work. 

Support of the Swiss National Science Foundation (project 20-63459.00) (J. P.M.), the Canadian 
Natural Sciences and Engineering Research Council (G.A.H.W.) and the National Research Council 
of Canada (D.A.B.) is gratefully acknowledged. 



- 7- 



Table 1. The Observations of ( Oph (HD 149757) 



3789 A 


5109 A 


K I Rad.Vel d 


T a S/N b W(1Q- 4 A) C 


T a S/N b W(W- 4 A) C 


(km s" 1 ) 


11700 2300 0.67 


10800 2600 0.50 


-14.53 ±0.18 



a Total exposure times in seconds. 
b per pixel 

c 3<r equivalent width detection limit for the band head (see section 2) 
d from Maier et al. (2001) 
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Fig. 1. — The final reduction sequence for the spectra at 3790 A. The flat field spectra contained a 
coarse structure absent from the stellar spectra which can be seen as common residuals in the dark 
subtracted and flat-fielded spectra for £ Oph and r\ UMa in the two bottom plots. The processed 
spectra of rj UMa were smoothed and then divided into the processed spectra of ( Oph thereby 
largely removing the flat-field structure. The final 3-point smoothed spectrum of £ Oph, the one 
discussed in the paper, is shown in the top plot. 
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Fig. 2. — Comparison of a laboratory spectrum (top) of C4 at 3788 A from Linnartz et al. (2000) 
smoothed to a spectral resolution of 110,000 and compared to the observed spectrum (lower) of ( 
Oph from Figure 1. 
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Wavelength (A) 



Fig. 3. — Comparison of a laboratory spectrum (top) of C5 at 5109 A from Motylewski et al. (1999) 
smoothed to a spectral resolution of 110,000 and compared to the observed spectrum (lower) of ( 
Oph. 



